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SUMMARY 


The problem of the oscillating lifting surface of finite span in . 
subsonic compressible flow is reduced to an integral equation. The 
kernel of the integral equation is approximated by a simpler expres- 
sion, on the basis of the assumption of sufficiently large aspect 7 
ratio. With this approximation the double integral occurring in the 
formulation of the problem is reduced to two single integrals, one of 
which is taken over the chord and the other over the span of the 
lifting surface. On the basis of this reduction the three-dimensional 
problem appears separated into two two-dimensional problems, one of 
them being effectively the problem of two-dimensional flow and the 
other being the problem of spanwise— circulation distribution. Earlier 
results concerning the oscillating lifting surface of finite span in- 
incompressible flow are contained in the present more general results. 

INTRODUCTION 


The present report is concerned with the problem of the oscillating 
airfoil of finite span, within the frame of the linearized lifting- 
surface theory. The aim of this study is the development of a theory 
which incorporates simultaneously the effects of three-dimensionality 
of the flow and of compressibility of the fluid. As an exact solution 
of this problem, even within the limitations of the linearized theory, 
presents very great difficulties, it is worth while to work toward an 
approximate theory which is valid provided the aspect ratio of the 
lifting surface is not too s m all. 

- The author has previously obtained results of this nature for the 
case of incompressible flow (references 1 and 2). In this earlier 
work the known results for the problem of two-dimensional incompressible 
flow were contained as a special case. The present work generalizes 
these results so as to take account of compressibility in the subsonic 
range. Thus the results of this report consist of a system of equations 
•which contain as special cases both the author* s results for the wing 
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of finite span in incompressible flow and the results of Possio’s theory 
of two-dimensional compressible flow (reference 3) . 

The scope of the present results may briefly be described as 
follows. The starting point of the work is an integral-equation 
formulation of the problem of the lifting surface of finite span. 

The integrals which occur are double integrals and the functions to 
be determined are functions of two independent variables. The essential 
step of the present work is to replace the actual kernel of the integral 
equation by an approximate kernel in such a way that the double integrals 
are reduced to single integrals over the range of either one of the two 
independent variables. In this way the problem is reduced to two 
problems which are to be solved separately. The first of these two 
problems is of the same nature as the Possio problem of two-dimensional 
compressible flow. The second of these problems is of the same nature 
as the problem of the Prandtl lifting-line theory for the wing of 
finite span in uniform motion. 

As in the theory of incompressible flow, this reduction of the 
double— integral problem to two single— integral problems depends 
crucially on the assumption of sufficiently large aspect ratio. While 
“sufficiently large" aspect ratios might be thought to be aspect ratios 
of about three, definite statements of this nature must be based on 
experimental evidence, as long as no exact solutions exist for the 
three-dimensional problem of the oscillating lifting surface in com- 
pressible flow. 

It is perhaps worth while to state explicitly that the present 
problem is quite, different from the corresponding problem for supersonic 
flow. 


It might also be added that there are reasons to believe that it 
is not satisfactory, even approximately, to superimpose aspect— ratio 
corrections for incompressible flow and compressibility corrections 
for two-dimensional flow in order to obtain corrections for the combined 
effect. This latter point is one of the reasons for the present study. 

In this first report the work is carried to the point where the 
double— integral equation is reduced to an equation containing single 
integrals only over the quantities to be determined. Further develop- 
ments will be given in a subsequent report. 

This work was conducted at the Massachusetts Institute of Technology 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 
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SYMBOLS 

Cartesian coordinates 

main— stream velocity in X— direction 

time 

defined by equation of lifting surface Z = H(X,Y,t) 

components of velocity change caused by presence 
of lifting surface 

region in X,Y— plane occupied by projection of 
lifting surface 

density of stream flowing with velocity U 

density and pressure changes, respectively, 

associated with velocity changes u, v, and w 

velocity of sound in main stream (a^ = dp/dp 0 ) 

coordinate of trailing edge of R_ 

Cl 

potential of velocity changes u, v, and w 
circular frequency of oscillation 
Mach number of main stream (M = U/a) 
real part of 

a length to be identified with the semichord 
of R a at midspan 

reduced— frequency parameter (k = cbb/U) 
dimensionless coordinates defined by equation (19) 
function defined by equation (21) 

kM^ 

parameter defined as p. = yz 

1 - yr 
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K 

parameter defined as k = — — 

1 - M 2 

V 

*|f 

parameter defined as v = 1 

1 -M 2 

V 

region in x,y— plane corresponding to region E a 
in X,Y— plane 


coordinate of trailing edge of E fl * 

X L 

coordinate of leading edge . of E a * 

V 

region in x,y— plane consisting of the strip to the 

right of the trailing edge of E fl * ^ 

E r * 

entire x,y— plane except for regions B a * and By* 

A 

function defined by equations ( 36 ) 

6 

function defined by equation (37) 

l,n 

variables of integration in accordance with 
equation (40) 

r 

defined by equation (4l) 


Hankel functions of second kind, and of zeroth and 
first order, respectively 

i 

auxiliary variable of integration 

K 

function defined by equation ( 51 ) 

G 

function defined by equation ( 52 ) 

In 

fimctions defined by equations ( 5 * 0 * ( 55 )* and ( 56 ) 

Pm 

order of magnitude of 

X 

function defined by equation ( 38 ) 

% 

function defined by equations (73) and ( 83 ) 
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auxiliary variable of integration 
auxiliary variable of integration 
function defined by equation (76) 
ratio of semispan to semichord at midspan 
defined as q* = tj/b \/ 1 — 

local semichord divided by semichord b at midspan 

a quantity indicating amount of sweep and defined by 
equation (86) 

dimensionless coordinate defined as x* = (x — xm*)/b* 
quantities defined by equation (87) 

quantities defined as g*(x*,y*) = g(x,y); X*(x*,y*) = X(x,y) 
function defined by equation (89) 
function defined by equation (90) 

THE BOUHDAKY— VALUE PROBLEM CSF THE 
OSCULATING LIFTING SURFACE 


It is assumed that a nearly plane, impenetrable surface is put into 
the path of an inviscid flowing fluid which, except for the effect of 
this surface, possesses a uniform velocity U in the direction of the 
positive X-axis. The impenetrable surface, henceforth called lifting 
surface, is taken to lie nearly in the X,Y— plane and its equation is 
written in the form Z = H(X,Y,t). When H = 0 no, disturbance is 
caused. When the lifting surface is not exactly plane and parallel to 
the direction of U the velocity components (U,0,0) are changed 
into (U + u,v,w) where u, v, and w depend on the form of the 
function H and on the shape of the region R a which is the projection 
of the lifting surface onto the X,Y— plane. 

The disturbances caused by the presence of the lifting surface 
are assumed to be small, in the sense that the differential equations 
and boundary conditions of the problem are linearized with respect to 
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the disturbance velocity components u, v, and w and with respect to 
the pressure and density changes p and p caused by the presence of , 
the lifting surface. 


Under these conditions the differential, equations of the problem 
are the following: 


c>U; + U — = d /P \ 
bt bX bX\po) 


( 1 ) 


+ u = - L7.E.N 

bt bX dYVpo/ 


, II _ _ d_/P\ 
bt bX bZ\P 0 J 


|e + tr |e + p y|“ ; * |tS . 0 

bt bX °\bX bY bzj 


( 2 ) 


(3) 

w 


p = a2p 


(5) 


The quant ity p Q in equations (l) to (4) is the density in the fluid 

flowing without disturbance and the quantity a in equation ( 5 ) is the 
velocity of sound in the undisturbed fluid, that is, a2 = dp/dp Q . 

The boundary condition of no relative normal flow at the lifting 
surface is satisfied, within the frame of the linearized theory, instead 
of on the lifting surface itself, on the projection of this surface onto 
the X,Y— plane. 


X,Y inside E a , w = + U g (6) 

The form of condition ( 6 ) (which holds on both sides of the lifting 
surface) indicates that w is an even function of Z. From equation ( 3 ), 
it follows then that p is an odd function of Z and the condition 
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that the pressure disturbance p is continuous, except when passing 
across the lifting surface^ means that for Z = 0 and 

X,Y outside R a , p = 0 ( 7 ) 


On the basis of conditions ( 7 ) and (6) the problem may be considered 
to consist in the determination of u, v, w, and p in the half 
space Z ^ 0 with the boundary Z = 0 . 

In addition to conditions (6) and ( 7 ) the following further 
conditions are prescribed in order to obtain an unambiguous solution. 

At the trailing edge of the lifting surface. 


X = It(Y), p is finite (8) 

Finally, it is postulated as a condition "at infinity" that energy is 
traveling outward without reflection, in a manner to be defined more 
precisely in what follows for the case of simple harmonic motion. 


VELOCITY— POTENTIAL FORMULATION OF THE PROBLEM 


It can be shown that the problem may be solved as stated in 
equations (l) to (8) by means of a velocity potential 0, in terms 
of which 


50 

U = — !■ 

5X 

50 

v 5 y 

w = M 

5z 


(9) 


Combination of equations ( 9 ), (l), (2), and ( 3 ) results in the following 
expression for the pressure change p: . 



( 10 ) 
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Combination of equations (10), (5), ( 9 ), and (4) results in the 
following differential equation for 0: 



Equation (ll) is to be solved in the half space Z > 0 subject to the 
following boundary conditions at Z = 0: 


X,Y inside. R a , ^ + u ^ 

a az at ax 


X,Y outside L. M + u M = 0 

a at ax 


X = Xjp(Y) , + U ^ is finite ( lit-) 


and subject to the condition of no energy .reflection "at infinity." 

In what follows attention is restricted to the case of simple 
harmonic motion and the following equation is written: 


( 12 ) 

(13) 


0(X,Y,Z,t) = e^CX^Y/z) 
with corresponding, expressions for H enrl p.^ 


(1“5) 


1 It is perhaps not entirely superfluous to indicate that this is 
meant in the sense that, corresponding to a surface equation ReCe^S), 
there is a pressure distribution ReCe^^p). 
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Equation (ll) now assumes the form 




bT 2 bZ 2 



+ 



= 0 


Equation (10) becomes 



(16) 


( 17 ) 


How the following dimensionless variables anti parameters are 
introduced: 





( 18 ) 


( 19 ) 


The differential equation (16) then assumes the form 
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Equation (20) is reduced further, for the purpose of eliminating 
f irsi>-derivative terms, by the following substitution: 


0 = e 1 ^ 


( 21 ) 


Choose 


d 



( 22 ) 


and obtain as the equation for i|r. 


+ K'njf = 0 


(23) 


where 



(24) 


How p and the boundary conditions must be expressed in terms of 
the new independent variables x, y, and z and in terms of the new 
dependent variable i|r. 

Prom equations ( 17 ) and ( 19 ) , it follows first that 

1 

< 25 

Combination of equations ( 25 ) and (21) gives 
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where V is defined as 



The boundary condition (12) becomes 


x,y inside Rg*, 


= E e itxx (ikH + |3) 

dz _ m2 \ Sx/ 


(27) 


(28) 


where R & * is the region in the x,y-plane corresponding to R a in 
the X,Yr-plane. 

The boundary condition (13) becomes 


x,y outside Ra*, 


ivty + ^ 
ox 


0 


(29) 


The condition of finite pressure at the trailing edge is now. 


x = s»r(y)^ 


ivi|r + ^ is finite 
dx 


(30) 


Finally, the condition of no energy reflection at infinity is written 
in the form 


z _^oo, \|r £ f(x,y,z)e i|tr (31) 

where r^ = x^ + y^ + z^ and where f tends to zero as z tends to 
infinity.^ 


^his ensures that 0 ~ fe i ( cofc-Kr ) as z -tends to infinity 
and therewith that waves are traveling away from the source of the 
disturbance . 
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The boundary condition (29) is made more specific in the following 
manner. From condition (29), 


*U,y,0) = c(y)e iVx (32) 

For any line y = Constant which does not pass through B a *, it can be 
concluded from the condition of undisturbed flow at X = —00 
that c(y) = 0. The same can be said for that portion of any 
line y = Constant which is situated in front of the leading edge 
of the airfoil region. The situation is different for portions of 
lines y = Constant which do pass through Eg* and which are to the 

rear of the trailing edge of E a *. The region to the rear of the 

trailing edge and bounded by lines y = Constant which are tangent 
to E a * is called the wake region and is designated by E^*. The 

exterior of the region R a * and E^* is called the remaining region 
and is designated by R r *. Then c(y) = 0 in Bj* and, in general, 
c(y) ^ 0 in Ry*. 

Equation (29) is thus seen to be equivalent to the following two 
equations: 


x,y inside R r *, t = 0 (33a) 

x,y inside B**, t = c(y)e _ivx (33b) 

There is now introduced a function A defined by 


A(y) = 2*[x r (y),y,C>] 


(3*0 


Equation (33b) can then be written in the form 


A(y)e iv ( x ^' x )’ 

dx 2 


x,y inside Rw*, 


(35) 
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In view of equations (33 a )j there may also he written 


A(y) 



dt(x,y, 0 ) 

dx 


dx 


( 36 ) 


where x^(y) is the coordinate of the leading edge of the airfoil 
region R a *. 


SUMMARY OF THE KESULTAHT BOUHBAHY— VALUE PROBLEM 


' Before proceeding with the solution of the problem as reduced in 
the foregoing section of this report its final formulation is reca- 
pitulated as follows. 

Determine the solution of the differential equation 

V 2 \|r + k 2 ^ = 0 (23) 

in the half space z > 0 subject to the following conditions at z = 0 : 


x,y inside R a *, = g(x,y) 

( 28 a) 

x,y inside R**, A(y) e iV ^^ 

(35) 

x,j inside R*,*, ^ = 0 

ox 

( 32 b) 

x = Xm(y), ^ is finite 
ox 

( 29 a) 


and subject to the following condition at infinity: 

z->», y = fe -iKr (31a) 

where f = 0(z) and where r^ = x^ + y^ + z^. 
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Various quantities occurring in these equations are defined as 
follows: 


g(x,y) 



A(y) 



X(x,y) dx 




2 dt(x,y, 0 ) 
dx 


(37) 
( 36 a) 

(38) 


The parameters k , M-* and v are defined hy equations (24), (22), 

and ( 27 ), respectively. The region R a * follows from the airfoil 

region R a hy multiplication with a scale factor l/b in the 

x— direction and hy multiplication with a scale factor Vl - M 2 /b in 
the y-direction. The region Ry* is the strip extending from the 

trailing edge x'= x/p to x = 00 , and the region Rj* is the 

remainder of the x,y— plane. 

The solution of the boundary— value problem is to be used to 
calculate the pressure— change amplitude p at the lifting surface 

in accordance with the relation 



(39) 


which follows from equation ( 26 ). 

The solution of the problem as summarized will be approached 
through its reduction to an integral equation for the quantity X as 
defined by equation ( 38 ). 
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AN INTEGRAL REPRESENTATION FOR THE VALUES OF dV/dz 


In this section it is proposed to derive a formula for the values 
of ch|r/dz, for the purpose of setting up the basic integral equation 
of the problem under consideration. To begin, results are taken which 
in essence are known and then are transformed in a way designed to 
facilitate the subsequent transition from the exact , double— integral- 
equation of the problem to the approximate integral equation containing 
single integrals only. 

The first formula is a representation of the values of ' cM/chc in 
the interior of the half space z > 0 in terms of the values of dt/dx 
on the boundary z = 0 of the half space, as follows: 

(40) 


In equation (kO) and in all that follows. 


r 2 = (x - |) 2 + (y - i)) 2 + z 2 (kl) 


and the quantity \ is according to equation ( 38 ) given by 

£,t), 0) /dZ. It is noted that this representation of dty/dx 
ensures that the conditions at infinity as expressed in equations ( 31 ) 
are satisfied. 

From equation (kO), it follows that 



If it is now observed that the quantity e — is a solution of 

the differential equation ( 23 ) for \|r, then equation (k2) may be 
written in the alternate form 


_ ^4 
dxdz k« 




(k3) 
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In equation (43), it is noted that = bP/br^ and the term in 

question is integrated "by parts with respect to tj. In addition to 
this, hy making use of the obvious identity 

X( t>Tj) =x(i,y) + [x(£,tj) -\(£,y)] 
equation (43) is written aB follows: . 



In the first integral on the right of equation (44) the integration 
with respect to t) may be carried out explicitly. When the remaining 
two integrals are absent, there is thereby obtained the appropriate form 
of the integral relation (i*) which would follow if two-dimensional flow 
had been assumed from the beginning. 

The following formula expresses the integral in question in terms 
of a Hanhel function 


—1 k\J( x-I) 2 +(y— tj) 2 +z 2 

O 

V(x - 0 2 + (y - tj ) 2 + z 2 


drj = iitH, 


: 0 ( 2 ) [kVu- O 2 + z 2 ] m 


The next step is 'to reduce the second integral on the right of 
equation (44) to an integral involving dX/dq. To this end, there is 
written 



- X(!,y)] 9 dtj 



U 


y 


(46) 
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The two integrals on the right of equation (1+6) are integrated "by parts 
and the constants of integration are chosen in such a way that the 
integrated portions vanish. After some elementary transformations, 
this leads to the following formula: 




U 



>-|) 2 +(y-Ti) 2 +z 2 
dn = 

/» 1 

\/( x - |) 2 + (y - Ti) 2 + Z 2 

^co 

|y - nl 

r i y e -tK\/(x-i) 2 +z 2 + 5 2 dt . 

5ii y-i 

J-to 

J - 00 - ^) 2 + Z 2 + t 2 


cLT) (1+7) 


Equations (1+5) and (1+7) are introduced into equation (41+) and the 
following relation is obtained: 


dxdz 


= J X(|,y)^^ + k 2 ^ j * H 0 (2 ) (k\|(x - I) 2 + z 2 ) d| + 

I y ^ -iK\/(x-|) 2 +Z 2 +^ 2 


1_ 

4n 




i y ~ ^1 /^- + K 2 ^ 


y - ^ 


V(x - I) 2 + z 2 + £ 2 


dt, -+ 


— iK\/(x— i)^+z 2 +(y— tj) 2 


^ V(x - 6) 2 + z 2 + (y - tj) j 




d| dq 


( 48 ) 


The fi nal step now consists in integrating equation (48) with 
respect to x between the limits — <» and x. The condition of 
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undisturbed flow far in front of the lifting surface makes (cM/8z) x _ _ 00 = 0. 
There is then obtained the relation 



Equation (49) expresses the values of dty /bz for z > 0 in terms 
of the values of X(g,q) = 2d'K|,rbO) fit and is thus the result, the 
attainment of which was the aim in the present section. 


THE INTEGRAL EQUATION OF THE OSCILLATING LIFTING SURFACE 


The integral equation of the problem is obtained by substituting 
the information contained in equations (28), (32), and (35) in 
equation (4-9) and by then letting z tend to zero. In so doing the 
integrals which contain infinities in the integrand are to be 
interpreted as Cauchy principal values. This latter step has been 
discussed in detail for the case of incompressible flow in reference 1. 
No additional difficulties in this regard are encountered for the 
present problem of compressible subsonic flow and the explicit 
justification for this step is therefore omitted in the present report. 
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The integral equation of the oscillating lifting surface in 
subsonic compressible flow is thus of the following form: 


g(x,y) - 


t»r(y) 


,(y) 


Ml,y)K[(x - (,)> k] as + 


f* A<y)e lv *rW 

2rt 


n°° 


e iV ^ k[(x — I); kJ dl + 


JxjXy) 



~ G [(x - |), (y - T|); k] d5 dt) - 
dt| 


iv 

kit 


|-(Ae iVX r)e iV ^ g[(x - |), (y - T}); k] d| dt) 




(50) 


The kernels K and G are of the • following form: 


' K *s(¥ b o < 2> ( |i “ 




^ ^( 2 )(|x« - l|) dx* 


(51) 
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Equation (50) holds when x and y are inside R & * and is to he 

solved for X, in terms of A and g, where A and g are defined 
by equations (36) and (37 ), respectively. 

For the case of two-dimensional flow, dX/dt) = 0 and the last two 
integrals in equation (50) are absent. The remainder of the present 
work has as its object the derivation of a procedure to take account of 
these last two integrals in an approximate manner which permits the 
calculation of the effect of three-dimensionality of the flow in a 
simpler way than by actually solving the complete equation (50). In 
the derivation of this procedure the integral over the airfoil 
region R a * and the integral over the wake region R w * are treated 

separately. 

REDUCTION OF THE DOUBLE HfTEGRAL OVER THE 

. • v 

AIRFOIL REGION R & * 

Write 
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where 



When dealing with lifting surfaces of sufficiently elongated 
form, that is, with surfaces of sufficiently high aspect ratio, 
jy — ,tj| » |x — 1 1 over the major portion of the surface. When this 
latter inequality holds, the terms I n can, as will he seen, be 
written in such a way that a relatively simple dominant term can he 
separated in each of them. These dominant terms will he used for the 
approximation to be developed. There remains then the question 
concerning the validity of the approximation over that portion of R a * 

where the inequality j x - « |y - tj| does not hold. This question 

is answered as follows. It is assumed that dx/dt) varies suff iciently 
slowly with T) so that in this portion of the’ region R a * it is 

effectively constant. If this is the case, all that is necessary is to 
take account of the fact that both the kernel G and the approximation 
to G to be obtained are odd functions of y — i), so that in both 
cases the contribution to the value of the integral coming from this 
portion of R a * can be neglected. 

The aforementioned argument is also implicit in the earlier 
derivations for the corresponding problem for incompressible flow 
(references 1 and 2). There appears to be no reason to, believe that 
this particular argument should be less applicable to the problem of 
subsonic compressible flow than to the problem of incompressible flow. 
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Equation (5*0 is written in the form 

n 


ii = — 


I 7 71 1 e -iiWU-0 2 +£ 2 |- x _ | 

« ; (x - 1) 2 + s 2 |_^ x - 1) 2 + 


+ 1 k(x - i) 




and this implies the following order-of— magnitude relations: 


II « -(x - 0 


M.-ikIcI/i x , 


IjL « -(x - |) 


od*±!^ ♦ 


jy-mv Vk-^r 

Equation ( 55 ) is written in the form 


Io = K C 


K 2 



f 1 S V-1 k\/(x , *) 2 +C 2 

+ p_ — — _ : — _ dx” 


—00 

Lvl 

Jo / \Z(*") 2 + C 2 


dC 


(57) 


P — I V— 1 


Io = K £ 


Ml 


^tHo ( 2 ) (k|C|) 


’x-4 


3 -1 k\/(i»») 2 .^ 2 

V(x **) 2 + £ 2 


dx 


»» 


dS (58) 


In equation (58), note that, when k is not too gm»n and 
when |y - q| » |x - || , 
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Px-| 


U 


0 


-i K\/(x”) 2 +f ; 2 

8 

tfU") 2 + ? 2 


dx’ * 



U-5) 


(60) 


Neglect the term in equation (60) compared with the term in 
equation ( 59 )> then. 


P-Ml 


Io £ 


4r B o (2) ( lt Kl) ^ 


(61) 


Note that this approximation ceases to be correct as k becomes 
smaller and smaller. However, when this is the case the contribution 
to the total coming from I 2 becomes negligible because of the 

factor k 2 in front of the integral in I 2 . 


Finally, write for I 3 




e -lK\l(x u ) 2 +(j-t [ ) 2 

Vu **) 2 + (y - n ) 2 


dx** 


I 


3 = 


3 in 
2 





e -iK\/(x ,, ) 2 +(y-ii ) 2 
V(x ") 2 + (y - T ]) 2 


(62) 


Now, when k | y — rj | is sufficiently large. 



(63) 
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and 


P_ 

Bt] 


P x ^ -iKN^x” ) 2 +(y— n) 2 

e 

V(x *’) 2 + 17 - Tl ) 2 


dx*' = Pm 


(x - l)e _iK l y_1 l| 


|y - n\‘ 


(x - £)k -i 

Hiy-nl 


-iK|y-T]| 


(64) 


Neglecting the terms in equation (64) compared vith the terms in 
equation ( 63 ), 


H ~ fc[f «o (2) ( K iy - ’ll)] - ¥" i r^ ! - h i ( 2> (M7 - n) ( 6 5) 


Comparison of expressions (57), ( 6 l) , and ( 65 ) shows that the 
contribution of Ip to the total may also he neglected. Introducing 

then expressions (6l) and ( 65 ) into equation (53), the following 
approximation is obtained 


f, 


j 


r 

n 

r 


^-ly-Til 

t &d 8 . d ’i~ 

L 

* 

a 

[ Px |y - Til 

2 

in 

Ptj y - 1 ) 

n, 

L- 

2 

—00 
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K jf h x ( 2 ) ( k Ij - 11 ) 


dl dtj 


(66) 


Equation (66) contains the fundamental simplification of the kernel G 
in the region B a *. 
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Noting that the factor of 8x/8q in equation (66) does not 
depend on |, there is written further 


^(t)) 


J x l( v ) 


8x 

C>T) 


dl 


8_ 

dq 




X dl - 


djCrp 


( X)xT ^ + (x)^ (67) 


d^ 

'*L dT) 


Take for (X) x its value zero immediately in front of the leading 
Xj 

edge and for (X)^ its value — ivA which follows from equations ( 35 ) 
and (38). Therewith, and with equations (36), equation (67) becomes 


\xr(n) 




Si 


dl = e- 1VBr(,l) ^( Ae - 1VIT ) 
dq 


(68) 


By introducing equation (68) into equation (66), the following equation 
is obtained: 


no 


u 

R a* 


8 q 
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irtK tx (2) 




dq 


( 69 ) 


Equation (69) represents the final result of the present section. 


REDUCTION OF THE DOUBLE INTEGRAL OVER 
THE WAKE REGION R^* 


In reducing the double integral over the airfoil region R & * to 
a simpler approximate form, use has been made of the fact that the span 
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of R Q * is appreciably larger than the chord of R & *. Evidently this 
is not the situation for the wake region Ry* and thus additional 
considerations are necessary for the integral extended over Ry.. 


Proceed as follows. Write the last integral on the right of 
equation (50) in the form 





d| 


drj 


( 70 ) 


The second of the two inner integrals on the right of equation (70) 
may be treated exactly like the integral over R & *. The first of the 

integrals has the property, as will be seen, that it is dependent 
on x in a simple explicit form. 


Taking the second integral first, there results, in analogy to 
equation (56), 
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— 00 
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e iv ^ d£ dr) 


( 71 ) 


Jx 


The remaining integral with respect to t in expression (71) will be 
introduced in evaluated form in the final collected form of the results. 



NACA TN 1953 


27 


Next the remaining integral in equation ( 70 ) is transformed, as 
follows 






It may he noted that the function reduces to the function F 

first introduced by Cicala when M = 0, that is, when the fluid is 
incompressible (references 4 and 5)* Combination of equations (72) 
and ( 52 ) results \ in the following form of F^, 
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Now -equations (71) and. (72) are combined in accordance with 
equation (70) in order to obtain the following approximation to the 
double integral over IV*! 
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(74) 
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THE APPROXIMATE INTEGRAL EQUATION OF THE PROBLEM 


Equations (7*0 and. (69) are substituted into equation (50) and 
the following approximate integral equation is obtained: 


g(x,y) 
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^Ue iVXfr )F M dq (75) 


Equation (75) represents in preliminary form the result to be 
obtained in this report. This result is reduced to a somewhat simpler 
form as follows. 

Set 


Ae 


iv *r,c 


(76) 


and introduce a new dimensionless spanwise coordinate y* defined by 


y * = y 

s\|l - M 2 


( 77 ) 
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If b represents the semichord at midspan and sb the semispan, 
then it follows from the definitions of x and y in equation (19) 
that the coordinate y* assumes values in the interval (—1,1) only. 

Further, set 

x L (y) = x L *(y*), »p(y) = VKy*) (78) 

and note that x L *(0) = -1, ^(0) = 1. For a rectangular lifting 
surface, - xl* = xp* = 1 throughout. 

From ’equations (51) and (65), there follows for the first integral 
on the right of equation (75) , 
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(79) 


The second integral on the right of equation (75) becomes 
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The third integral, on the right of equation (75) becomes, 
with k from equation (24), 



n 
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(81) 


Finally, there is written for the fourth integral on the right of 
equation (75 ) j 


^-(Ae 1VXr )F M (y - tj) drj ■ 
dTj 
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dry* 
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(82) 
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where the function as defined by equation (73), may be written 

in the following alternate form: 


F„(x) - M 
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Now collect equations (82), (8l), (80), and (79) and substitute 
the result in equation (75) • This gives 


g(x,y*) = 
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Finally, there are introduced, dimensionless coordinates x* 
and I* as follows: 


x* = 


2x - (xr + x L ) 
*r - X L 


( 85 ) 


so that the interval x L < x < aj. goes over into the 
interval -1 ^ x* ^ 1. Further, 


Xiji — Xj^ = 2b* 

f PT + *l) = V 


( 86 ) 


where b* is the ratio of local ' semichord to semichord b at midspan. 
Finally, the following quantities are introduced: 


k* = b*K 


v* = b*v > 
ft* = ty/b* 


( 87 ) 


Note that, for the rectangular plan form, 
the elliptical plan form, b* = \jl — y*^. 


b* = 1, xjn* = 0 and, for 

V = °* 
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Equation (84) then assumes the following form: 
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Comparison of equations (88) and (84) gives for the kernels N 
and Sjj the following expressions: 
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The taslj from here on is the following. Equation (88) must "be 
solved for A*, in terms of g* and ft*. This part of the problem 
is exactly as in the two-dimensional theory. The function ft* is 
then to he determined by an integrodifferential equation which is 
obtained by expressing ft* in terms of A* in accordance with the 
definition of ft*. This part of the problem is similar to the earlier 
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work on incompressible flow of references 1, 2, and 6 . Finally the 
solution of the integral equation for ft* must be used to obtain 
expressions for the pressure distribution at the airfoil, as affected 
by the three-dimensionality of the flow about a wing of finite span. 

The results, as expressed by equations ( 88 ), ( 89 ), (90)* and ( 83 ) 
include the special case of two-dimensional flow for which dft*/dT)* = 
and the special case of incompressible flow for which M = 0. They 
also include essentially known results on compressible steady flow for 
which k = 0 . 


Massachusetts Institute of Technology 
Cambridge, Mass., May 18, 19^8 
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